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Abstract This work studies the concentration of radon in soil around a fault 
in the East Franklin Mountains in the El Paso area in West Texas. It is found 
that the in-soil production of radon is correlated to the existence of a fault 
even if it has not had any recorded activity in recent geological times. This 
adds to previous observations that link the production of radon to seismic 
activity, and seems to indicate that in non-active faults the radon production 
is due mainly to the radioactivity of the top soil and to the transport properties 
of the medium and not to deeper seismic activity. These results open the 
possibility of using in-soil radon gas concentrations as an examination tool 
of dormant faults. 
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1. INTRODUCTION 


Radon is a gas that is produced as part of the radioactive decay of 238 U, which 
includes 23 °. 234 Th, 234 Pa and 234 U and other elements that exist abundantly in 
most types of rocks and sediments; radon is thus present in all types of soil. 
The concentration of underground radon near the surface depends solely on 
the presence of its parent nucleus, 226 Ra, and the transport properties of the 
medium. As the medium transport properties vary drastically near geological 
faults, it is expected that in-soil radon gas concentration will increase near 
faults [5]. 

The increase of radon concentrations in the neighborhood of active faults, 
has been explained in terms of the dilatancy-diffusion model of earthquake 
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occurrence [21], and has even been used to predict tectonic activity. In a 
nutshell, changes of elastic strain before earthquakes cause rocks to dilate 
forming open fractures and facilitating rupture, which in turn allows the flux 
of gases including radon. Since radon has a very poor intrinsic mobility in 
soil and a very short half-life (3.8 d), it would not be expected to propagate 
upwards long distances, however, the enhanced upward speed of the radon 
gas is explained by the existence of rising fluids such as underground water, 
carbon dioxide, nitrogen and methane produced and liberated by the tectonic 
activity [4]. Such condition of increased flow is expected to exist for as long 
the fault continues being active. 

Due to wide variation of conditions, no precise world averages exist for 
in-soil concentrations of radon nor for its variation around fractures and faults. 
Increments ranging up to 3.233% over the mean background radon level have 
been reported at the Bad Nauheim fault in Germany [12], where a peak radon 
signal of over 1,000 kBq/m 3 was observed around a region where the mean 
background level is 30 kBq/m 3 . On the other hand, since similar increments of 
radon concentration have also been reported around less active faults [4, 14, 
22], it is possible to think that enhanced radon exhalations can occur in even in 
faults without any recent tectonic activity. 


2. HYPOTHESIS 

In this article we test a complementary hypothesis, namely, that once a fracture 
or fault has been created, the enlarged path for gas flow remains and facilitates 
the flow of radon gas even in the absence of any tectonic activity. That is, 
in-soil radon concentrations should increase near faults as compared to the 
immediate background concentration of the neighborhood. 

It is clear that, without the help of tectonic gases, the underground radon 
gas found near the surface could not have been produced at greater depths. 
Radon, at a difference from tectonic gases, is produced at all depths, but due to 
its slow diffusion and short half-life, it can only migrate short distances. 

To understand this it is necessary to remember that 222 Rn is generated 
through the 238 U decay series which produces several long half-life isotopes 
( 234 U, 230 Th and 226 Ra) which are commonly found in granite, igneous rocks, 
sedimentary rocks, metamorphic rocks, etc. It is this abundance that permits 
the production of 222 Rn by practically all types of rocks and soils at all depths. 

Furthermore, the radon ascent velocity in soils is relatively slow, and its 
mean distance displacement is relatively short due to its brief half-life. For 
example, taking into account the 222 Rn half-life of 3.82 d and the diffusion 
coefficient of radon in dry soil (5x 10' 6 m 2 /s), the average diffusion length is 
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found to be 1.6 m [15]. Without tectonic gases, radon gas produced at much 
larger depths would not live long enough as to reach the surface through dry 
soil. The existence of fractures and faults, however, could increase the diffusion 
length of gas produced near the surface allowing a much larger concentration 
of gas around faults. 

As a proof of principle, we test whether radon concentration increases 
noticeable near dormant faults, i.e. around faults that have not had any recorded 
activity in recent geological times. In this article we present a preliminary 
study of 222 Rn measurements along a dormant fault in the Franklin Mountains 
in the El Paso Area in West Texas, and contrast it with studies that identify 
fractures and faults through the general area. If proven correct, this hypothesis 
would open the possibility of using in-soil radon gas concentrations as an 
examination tool of dormant faults. 
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3. GEOLOGICAL SETTING AND LAULT TRACE LOCATION 

Measurements of radon concentrations were taken on the east side of the 
Franklin Mountain in El Paso, Texas in the United States of America. These 
are tilted block-fault mountains formed by fractures of the Earth’s crust 
produced by near-vertical faults which thrusted and tilted the landmass upward 
producing diagonal rock layers which expose layers of Precambrian rock 1.25 
billion years old. The area of study (see Figure 1) is located at the east flank 
of the Franklin Mountains, which constitutes a small range trending north- 
south between south New Mexico and the western tip of Texas. Although the 
Franklin Mountains are mainly composed of Cretaceous sedimentary rock 
units formed when the region was covered by water, there are also tertiary 
igneous intrusions and Precambrian crystalline rocks. 

Tectonically, two mayor pulses, the Laramide orogeny (85-45 My) and the 
Cenozoic Rio Grande Rift, have shaped the geological landscape around the 
area of study. Flowever, the major structural control on the uplift and tilting of 
this range is the Late Cenozoic crustal extension process associated with the 
Rio Grande Rift, which has resulted on a well-defined horst and graben classic 
rift architecture. The Bolson del Flueco, located at the eastern upskirts of the 
Franklin Mountains is an extensional basin developed as consequence of the 
Neotectonic displacement of the basement by more than 2,700 m in depth 
along a series of quaternary age of the Late Miocene [13]; there are normal 
faults bounding the basin between the Franklin Mountains and the Hueco 
Tanks [6]. 

The main fault evidenced by surface geology is the East Franklin Mountain 
Fault (EFMF), which has been dated as 8-12 ky [16], with scarps ranging 
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Figure 1: Location of the Franklin Mountains, the area of study, collected points 
and approximate location of the East Franklin Mountain Fault (EFMF). The inset (a) 
shows the tilt derivative map derived from the Corrected Bouguer Anomaly to infer 
the trace of the EFMF. 


between 2 m and 60 m [16, 6], and with average displacements of 0.18 mm/ 
yr [17]. This 52.7 km length fault [17] is considered, assuming and equivalent 
surface rupture, with potential to trigger a 7.1 Richter scale earthquake ([23], 
normal fault data). The EFMF is also clearly evidenced by the strong gradient 
on the tilt derivative map of the gravity corrected Bouguer anomaly (Figure 
1, inset a), derived from the gravity data base of the University of Texas at El 
Paso. 
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Figure 2: Fault trace location and 2D profiling. Modified from [3]. 


Although the East Franklin Mountain Fault surface location is mapped and 
available through the USGS quaternary fault database [13], the fault attributes 
are not detailed. Fault anatomy can be found on more recent studies which have 
not only mapped the quaternary faults in the El Paso Texas - Ciudad Juarez 
metroplex area, but also modeled gravity 2D profiles [2, 3, 18, 7]. The study of 
Avila [3] shows two profiles located nearby the study area of this research. Figure 
2, modified from such study, shows the fault trace location and 2D profiling 
of the zone. In particular profile B-B', located 6 km north to the radon field 
collection area, shows how the gravimetric signature can be properly fitted with 
a 2000 m heave, and a trow of nearly 2600 m. This quasivertical displacement of 
the Precambrian basement is associated with the EFMF. Avila et al. also shows 
the gravity forward modeling correlated to a north located seismic reflection 
section [8] providing both geophysical and geological constraints that support 
the normal fault controlled horst-and-graben basin architecture. The gravity 
forward model also reveals the outcropping Precambrian basement that might 
be acting as radioactive source. 

In this study we connected the EFMF along Avila’s profiles BB and CC. 
Since the radon decay products are expected to occur across faults and nearby 
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presence of crystalline rock units, the fault trace was mapped and identified 
on the terrain to collect a series of radon concentration readings on the field as 
shown in Figure 1. 

4. RADON MEASUREMENTS 

Although the natural radon family has among their members the isotopes 222, 
220 and 219, the most often measured is the radioactive inert gas 222 Rn that 
has a half-life of 3.82 days and belongs to the 238 U chain. The world average 
exhalation flux at the soil surface is 22 mBq/m 2 s which corresponds to about 
40 Bq/m 3 indoors; however depending on the subsurface gas permeability, the 
flux can be much larger. 

Measurements of radon content in ground air were carried out using a 
Markus 10 portable instrument manufactured by Gammadata in Sweden. The 
instrument operates by pumping out gas from the ground through a probe 
inserted in the soil at a relaxation depth of ~35cm during 30s. The Rn gas 
is stored for 10 min. in a measuring cell, time during which it decays into 
its progeny ( 218 Po and 214 Po) through alpha radiation detected with a silicon 
surface barrier detector operating under an strong electric field (~600 VDC). 
The readings are converted by the device into volume concentrations (i.e. an 
average radon concentration in the air-filled part of soil-pores) in kBq/m 3 . 
Table 1 shows the radon measurements obtained at the 16 points selected. 

The nominal measuring error established by the factory calibration is of 
1.7% of the measurement, but the readings are additionally rounded off by the 
LED readout of the device to tenths of kBq/m 3 ; this sets the average uncertainty 


Table 1. Location of points and radon concentrations obtained. 


Point 

Easting Northing 

(m) (m) 

Radon 

Concentration Point 
(kBq/m 3 ) 

Easting Northing 

(m) (m) 

Radon 

Concentration 

(kBq/m 3 ) 

A 

361731 

3525816 

1.05 ±0.05 

I 

361553 

3526166 

2.75 ± 0.05 

B 

361883 

3525823 

8.05 ±0.137 

J 

361555 

3526062 

3.75 ± 0.064 

C 

361586 

3525789 

8.75 ±0.149 

K 

361606 

3526056 

1.35 ±0.05 

D 

361660 

3525801 

4.05 ± 0.069 

L 

361819 

3526143 

4.35 ± 0.074 

E 

361436 

3525831 

4.35 ± 0.074 

M 

361706 

3526041 

1.75 ±0.05 

F 

361286 

3525855 

5.75 ± 0.098 

N 

361787 

3526655 

7.75 ±0.132 

G 

361716 

3526156 

5.75 ± 0.098 

O 

361830 

3526632 

6.05 ±0.103 

H 

361632 

3526159 

9.35 ±0.159 

P 

361733 

3526659 

7.05 ±0.120 
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to a minimum of 50 Bq/m 3 and as high as 1.7% of the measurement. For this 
reason we take all readings as being in the middle of each 100 Bq/m 3 bin and 
list the uncertainty as the maximum of 1.7% or 50 Bq/m 3 . For instance, the 
reading at point C, which was recorded by the Makus 10 as 8.7 kBq/m 3 , could 
indeed correspond to any value between 8.700 and 8.799 kBq/m 3 , thus in Table 
1 it is listed as 8.75 kBq/m 3 with a 1.7% error of 0.149 kBq/m 3 . 

The first observation is that the overall scale of measurements is in the 
same range as observed in many active faults as reported, for instance, by 
[19] (see Table 1); indeed these readings are higher than those obtained near 
the active fault of Dead Sea Transform in Jordan (1.8 - 3 kBq/m 3 ) [1] and 
the Johnson Valley Fault in Landers, California (4 kBq/m 3 ) [19], similar to 
those on the North Anatolian Fault in Turkey (9.8 kBq/m 3 ) [10], and not too 
distant to those measured on the North and Northwestern Greece Fault (13 
kBq/m 3 ) [11]. This immediately appears to support the hypothesis that radon 
production is not related to seismic activity. 

Overlaying the values of concentration of Table 1 over the geographic 
map of the EFMF it is possible to obtain a contour map of the concentration 
levels of the inspected area; Figure 3 shows such contour map using a color 
scale ranging from 1.1 kBq/m 3 to 9.1 kBq/m 3 . Although the lack of information 
at the far edges of the map distorts a bit the information presented, the map 
clearly shows that the highest concentration levels, namely points C (8.7 kBq/ 
m 3 ), H (9.3 kBq/m 3 ) and N (7.7 kBq/m 3 ), coincide with the expected location 
of the EFMF Furthermore, the concentrations are observed to die down very 
rapidly for points merely meters away from the fault. 

The noticeable exceptions are the points J, K and M which show much 
lower levels of radon concentration: 3.7, 1.3 and 1.7 kBq/m3 respectively. 
These points, fortuitously, are located on a dirt road that was excavated 
next to a wash which, undoubtedly, removed and eroded soil and deposited 
extraneous debris of several feet of depth, that affected both the radon 
production (compared to neighboring terrains) and the structure of the 
uppermost soil and its gas transport properties. Seen from the positive side, 
the low readings obtained at the dirt road and wash seem to confirm that the 
in-soil radon concentration is produced mostly by the top soil and not by the 
deeper undergrounds which, in this case, were covered by dirt and debris of 
the erosion. 

[A second possible conclusion, not explored in this article, is that the 
reduction of radon gas could indicate a segmentation of the fault. This would 
be welcome news for the El Paso Area population as a 52.7 km fault carries a 
risk of 7 in the Richter scale while two faults of 26 km would reduce the risk 
to 5 in the same scale; further analysis would be needed to test this possibility.] 
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STUDY AREA LOCATION 
Radon measures 
points 


Quaternary faults 

_(United States 

Geological Survey) 
Radon measures 


Figure 3: Concentration of radon gas measured on the area studied; the color scale 
varies from 1.1 kBq/m 3 (blue) to 9.1 kBq/m 3 (pink). The solid line approximates the 
location of the EFMF. 


As pointed in other studies [9], the soil-gas concentrations of radon are 
expected to vary in time due to several factors (temperature, precipitation, 
humidity, pressure, tectonic activity), thus, to fully substantiate the hypothesis 
presented in this article a continuous long time monitoring of this and other 
faults is necessary. It must be pointed out that, to support the hypothesis, the 
levels at the fault would only have to peak with respect to the background of the 
immediate neighborhood and not attain specific levels. Such spatiotemporal 
study is currently underway [20]. 

5. CONCLUSIONS 

By studying the production of radon in soil around the East Franklin Mountain 
Fault in El Paso, Texas, we determine a strong correlation between the in¬ 
soil concentration of radon gas and the location of the fault. The fact that 
the fault has not had any activity in recent geological times serves as a proof 
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of principle that enhancement of radon production near a fault is not due to 
seismic activity as predicted by the dilatancy-diffusion model of earthquake 
occurrence, but to the inherent radioactivity of the soil and to the enhanced 
gaseous transport properties of the soil due to the existence of the fault. Since 
the in-soil radon concentration is observed to decrease where extraneous dirt 
and debris have been deposited on top of the fault, this supports the idea that 
the in-soil radon concentration is produced mostly near the top soil and not 
deeper in the ground. 

We conclude that the use of non-invasive measurements of radon gas 
concentrations to detect non-active faults appears as a real possibility that 
deserves further study. If proven correct by further studies, this method of fault 
detection has distinct advantages when used in populated areas. For instance, it 
does not rely on the production of shear and compressional waves, such as the 
method of seismic tomography, nor needs to produce electrical signals which 
can interfere with the urban electrical network as does the method of electrical 
resistivity tomography. These advantages are especially useful in determining 
the location of faults inside cities. 
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